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By considering three different spin models belonging to the generalized voter class for ordering 
dynamics in two dimensions [I. Dornic, et al. Phys. Rev. Lett. 87, 045701 (2001)], we show that 
they behave differently from the linear voter model when the initial configuration is an unbalanced 
mixture up and down spins. In particular we show that for nonlinear voter models the exit proba- 
bility (probability to end with all spins up when starting with an initial fraction x of them) assumes 
a nontrivial shape. The change is traced back to the strong nonconservation of the average mag- 
netization during the early stages of dynamics. Also the time needed to reach the final consensus 
state Tn(x) has an anomalous nonuniversal dependence on x. 

PACS numbers: 05.40.-a, 89.65.-s, 64.60.De 



I. INTRODUCTION 

The voter model (VM) [TJ [5] is a paradigm of coarsen- 
ing phenomena [3] that stands as one of the most interest- 
ing models in non-equilibrium statistical mechanics 4J. 
The nature of its interest is twofold: On the one hand, 
it represents one of the few nontrivial non-equilibrium 
statistical processes that can be exactly solved in any 
number of dimensions On the other hand, it has a 

natural application in social dynamics [7J as a model for 
the formation of opinion consensus in a society initially 
divided in two different standpoints. The appeal of the 
VM is further enhanced by its connection with neutral 
models in genetics, ecology and linguistics [SHTU], Its def- 
inition is very simple: On a regular lattice or graph, each 
site is endowed with a binary variable = ±1. At each 
time step, a randomly chosen site copies the state of one 
of its nearest neighbors, chosen in its turn at random. 
This parameter-free dynamics can be succinctly encoded 
in the flipping probability f(xi), measuring the probabil- 
ity that spin i will flip if surrounded by a fraction Xi of 
spins in the opposite state, which takes the simple linear 
form f(xi) = Xi. For this reason we will refer to it in 
the following as the linear voter model. Voter dynamics 
is thus characterized by the presence of two absorbing 
states (all spins either +1 or —1, the consensus states) 
with a Z2 spin reversal symmetry. Moreover, since the 
rate of creation of +1 and — 1 spins is equal, the magne- 
tization is conserved in average. 

The way in which consensus is reached in the VM can 
be characterized from different perspectives. From the 
point of view of non-equilibrium statistical mechanics, 
the coarsening process in the VM is marked by the ab- 
sence of surface tension causing an anomalous loga- 
rithmic decay of the density of interfaces in d = 2, namely 
p(t) ~ 1/ ln(i), in opposition to curvature-driven dynam- 
ics [3], which leads to an algebraic decay pit) ~ V 1 ! 2 . In 
the social dynamics context, on the other hand, interest 



is focused on the exit probability E(x) [7J (defined as the 
probability that the final state corresponds to all sites 
in state +1) and the consensus time Tjf(x) (the average 
time needed to reach consensus in a system of size N) 
when starting from fully random initial conditions with 
a fraction x of sites in state +1. Conservation of mag- 
netization implies a characteristic linear exit probability, 
E(x) = x, in any dimension al [3], while the consensus 
time takes the form, for d > 1, [T^] 

T N (x) = -N cS [x\n(x) + (1 - x)\n(l - x)}, (1) 

where N e g is an effective factor depending on the number 
of sites N and dimensionality d [I]. 

A detailed analysis revealed that VM actually lies at 
the transition point between a ferromagnetic (ordered) 
and a paramagnetic (disordered) phase, such that in- 
finitesimally small perturbations are able to drastically 
change its behavior [T3TU5| . The parameter-free nature 
of the VM thus led naturally to the question as whether 
it represents a peculiar and isolated point, or rather be- 
longs to a more general (universal) class of models, shar- 
ing the same properties. This issue has been answered 
by Dornic and coworkers [HI US] (see also [U1HE]), 
have pointed out the existence of a genuine generalized 
voter (GV) universality class, encompassing systems at 
an order-disorder transition driven only by interfacial 
noise, between two "dynamically equivalent" absorbing 
states. The dynamical equivalence between states can be 
enforced either by Z 2 -symmetric local rules, or by global 
conservation of the magnetization. Linear voter model 
possesses both properties, but each of them separately is 
sufficient to ensure GV behavior. The GV class is char- 
acterized in d — 2 by the logarithmic decay of the density 
of interfaces F] as well as by other critical exponents [TT] . 



The behavior of the VM in d = 1 coincides with the zero temper- 
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In this generalized perspective, the GV transition for Z 2 - 
symmetric models can be theoretically rationalized as the 
superposition of two independent transitions [TBI HU > an 
Ising transition and a directed percolation |20j transi- 
tion, whose respective symmetries are broken in unison 
at the GV manifold. These ^-symmetric models are also 
called nonlinear voter models, because at the transition 
the flipping probability f{xi) assumes a nonlinear form. 

In this paper we present a detailed numerical investiga- 
tion of three types of models belonging to the GV class: a 
non- linear voter (NLV), originally devised to explore the 
GV manifold the recently proposed non- linear q- 

voter model (qV) (2JJ; and the Kaya, Kabakcioglu, and 
Erzan (KKE) model [22]. By means of extensive nu- 
merical simulations, we confirm the existence of the GV 
universality class, which is characterized in terms of the 
logarithmic decay of the density of interfaces at the GV 
point as well by the values of critical exponents related 
to the fluctuations of the magnetization and the corre- 
lation length when approaching the critical point from 
the disordered and ordered phases, respectively. In this 
respect, our simulations suggest non mean-field correla- 
tion exponents, at odds with the claim made in Ref. 
Surprisingly, when probing the behavior of the models 
for unbalanced initial conditions (x 7^ 1/2) by means of 
the dependence of the exit probability and of the con- 
sensus time on x, we observe that the originally defined 
GV class exhibits however strong nonuniversal features. 
In particular, we find that, at criticality, quantities such 
as the exit probability and the consensus time depend 
on further microscopic details of the models undergo- 
ing the GV transition. We show instead that models in 
which conservation of average magnetization is strictly 
enforced, such as the KKE model, have indeed a linear 
E{x) as the VM, but they exhibit a consensus time Tjv(x) 
different from the entropic form Eq. ([I]). On the other 
hand, models which exhibit ^-symmetry display E{x) 
and Tjy(x) both departing from the VM behavior. The 
non-linearity of the exit probability can be rationalized 
by inspecting the behavior of the average magnetization 
over time. Here we can see that magnetization is not con- 
served over short time scales, but it increases initially in 
the NLV model, while it decreases in the qV model. This 
transient behavior can be traced back to the presence of 
strong non-zero drift in the initial dynamical evolution 
starting from x 0. After this initial drift vanishes, the 
average magnetization remains constant and the ensuing 
evolution is well described by a linear voter dynamics. 

We have organized the paper as follows: In Section [IT] 
we present the results of numerical simulations for the 
class of non-linear voter models, determining the values 
of the critical exponents and showing the nontrivial x 



ature Glauber dynamics [4], while it is described by mean- field 
theory for d > 2, its upper critical dimension. The interest of its 
definition and properties is thus essentially given by the behavior 
in d = 2. 



dependence of the exit probability and of the consensus 
time. In Secti on|III| we do the same for the q- voter model, 
while Section IIVI is devoted to the KKE model. The 
final Section summarizes the results and discusses their 
relevance. 



II. INTRODUCING THE GENERALIZED 
VOTER CLASS: THE NON-LINEAR VOTER 
MODEL 

The characteristics of the GV class were exposed in 
Ref. by numerical examination of a Z 2 -symmetric 
non-linear voter model (NLV) defined it terms of a ki- 
netic Ising model as follows: We consider a binary spin 
system in d = 2, in which the probability r s ./j that a spin 
s flips depends on its value and the value of the local 
field h it feels. The Z2 symmetry imposes r_ s ,h = r Sj _^,; 
therefore all flipping probabilities can be encoded in the 
flipping probability for a s = +1 spin, r/,. = The 
absence of bulk noise imposes = 0. The standard 
linear VM is given by = 1/2 — h/8. In the general 
case, the NLV model depends on four free parameters, 
r_4, r_ 2 , ro and r 2 . In the following, we adopt the arbi- 
trary parametrization of Ref. |llj . imposing r_ 2 = 0.275, 
ro = 1/2, r 2 = r_4/4, and taking r_4 = e as a free tuning 
parameter. 

With this parametrization, the GV point corresponds 
to a critical value e c separating a paramagnetic phase for 
e > e c from a ferromagnetic phase at e < e c . The critical 
point can thus be estimated by monitoring the density 
of interfaces p(t) and identifying e c as the value leading 
to a logarithmic decay p(t) ~ l/ln(i), separating a con- 
stant behavior (e > e c ) from an algebraic decay (e < e c ) 
|11) . Here we follow a different approach, based on the 
exit probability E(x). The investigation of this quantity 
represents an effective way to determine the critical point 
e c . Indeed, e > e c corresponds to a disordered paramag- 
netic phase with E(x) — 1/2, while e < e c corresponds 
to an ordered ferromagnetic phase, where E(x) = Q(x), 
the Heaviside theta function. Therefore, focusing on an 
initial density x < 1/2, we should observe E(x) — > 1/2 
for e > s c , E(x) for e < e c , and E(x) — > const < 1/2 
for e = s c when increasing the system size L. 

In Fig. [TJ (main plot) we report the exit probability 
for x — 0.25 and different values of £ as a function of the 
lattice size L. A plateau is obtained for e ~ 0.3996, while 
larger (smaller) values of e lead to an increase (decrease) 
of E(x — 0.25) with L. We conclude that the critical 
point of the GV transition in d = 2 for the NLV model 
is located at e c = 0.3996(4), in good agreement with 
the result inferred from Fig. 3(a) in Ref. [11) . where a 
behavior for the density of interfaces with the logarithmic 
VM form was found. 

The properties of the GV universality class can be fur- 
ther explored by considering several critical exponents, 
measured in the vicinity of the critical point e c . These 
exponents are usually defined in terms of the suscepti- 
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FIG. 1. (Color online) Main: Finite-size scaling of the exit 
probability at x = 0.25 as a function of system size in the NLV 
model in d = 2, computed over at least 5 x 10 5 independent 
realizations of the dynamics for each system size. Inset: Sus- 
ceptibility and correlation length as a function of A = \e c — £ 
on lattices of size L = 100 and L = 5000, respectively. 



bility, measured as the fluctuations of the magnetization 
= J2t S i/N, i.e. 



(2) 



when approaching the transition from the paramagnetic, 
disordered phase, and the correlation length £ which, 
when approaching the GV manifold form the ferromag- 
netic, ordered phase, can be measured from the relation 

Em 

p(t)~#-V2. (3) 

Close to the critical point, these two quantities depend 
on A = \e c — e\, defining the critical exponents 

X (A) ~ A -7 , £(A) ~ A"". (4) 

In Fig. [I] (inset) we present the results of numerical 
simulations of the quantities £ and % as a function of 
A. While the determination of 7 is straightforward, the 
measurement of v is hindered by extremely long pre- 
asymptotic effects in the curvature-driven regime, leading 
to a decay of p(t) with an effective numerical exponent 
smaller than 1/2 [23] • Here, in order to obtain informa- 
tion about £, we proceed by performing a linear regres- 
sion of l/p 2 (t) as a function of t, and assigning to £ the 
value of the slope thus obtained. Data obtained in this 
way (Fig. [T] (inset)) provide the exponent values v ~ 0.60, 
7 ~ 1.26. The value of 7 is in excellent agreement with 
early numerical values |11( 1 1 3] , while v is rather different 
from the estimate of Dornic et al. [TT] and compatible 
with the result from [13] . Both exponents are also quite 
compatible with the scaling relation 7 = 2v. With re- 
spect to the mean-field values (v = 1/2, 7 = 1, with 
logarithmic corrections) proposed in Ref. |llj . from our 
data it is difficult to make a definite discrimination for 



the exponent v, since the plot of £ as a function of A 
can be equally well fitted to a pure power-law with non 
mean-field exponent or to a mean-field value with log- 
arithmic corrections. On the other hand, the exponent 
7 seems apparently better fitted with a non mean-field 
power-law exponent. 

To check these results we consider the linear VM in 
the GV manifold, which in the NLV model introduced in 
[TTj can be approached by setting r_4 = 1, ro = 1/2 and 
taking 7-_ 2 = e —> e c = 3/4. Here the state is paramag- 
netic for e < e c , while it is ferromagnetic for e > e c . In 
this case, see Fig. [I] (inset), we obtain 7vm — 1-29 and 
z^Vm — 0.62, which confirm the universality of the GV 
class. 

The analysis presented above confirms the results of 
previous studies. However, it also points out a new 
and surprising feature, which is only evident in sim- 
ulations performed out of the initial symmetric state 
(x = 1/2). As we can see from Fig. [I] (main plot), the 
exit probability of the NLV model at criticality (on the 
GV manifold) computed at the non-symmetric homoge- 
neous initial state x = 0.25, takes a value E(x = 0.25) = 
0.32 ± 0.02, i.e. larger than 0.25, well beyond error bars. 
This observation hints towards a non-linear form of the 
exit probability which is confirmed in Fig. pk). In this 
plot we can see that the exit probability E(x) deviates 
from linearity for the whole range of values of x. This 
deviation from linear VM behavior extends also to the 
consensus time as a function of x, as we can also see in 
Fig. §3). 

This departure of the NLV model from linear VM be- 
havior can be understood by looking in detail at the 
time evolution of the magnetization <f>(t) in the system, 
starting from x(0) — 0.25, corresponding to (f)(0) — 
2x(0) — 1 = —0.5, see Fig. [3j Data shows that mag- 
netization is strongly not conserved at short times, but 
in fact it experiences a sharp increase until it stabilizes, 
for times t > 100, at a plateau with approximate value 
4>cg ~ —0.37. The peculiar time evolution of the average 
magnetization (already noticed in Ref. |llj ) can be re- 
lated to the drift v((f>) in a Langevin representation [16j . 
in the form d t ((f>) = (v((f>)) [23]. We estimate the average 
drift (v((f>)) by computing — 2/N J2i s if( x i) (where Xi is 
the fraction of neighbors of i in opposite state) and we 
average all values of drift with the same magnetization 
value (f> = l/NJ2i s i- I n Fig. [4]we plot (v(<f>)) vs 4> for 
different system sizes in two distinct temporal regimes. 
For short times (t < 100) a sharp rise is present in the 
vicinity of the initial magnetization. This is responsible 
for the initial increase of magnetization until it reaches 
the steady state and its conserved (in average) value. For 
larger values of time t > 100 this rise is absent, and the 
drift takes a flat, almost vanishing form, thus ensuring 
conservation of magnetization. 

Given the flipping probability f(x), the origin of the 
rise for short times can be understood by considering the 
initial uncorrelated condition. In that case the average 
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FIG. 2. (Color online) (a) Exit probability as a function of 
the initial density for different system sizes in the NLV model 
and the q — 4 qV model at the respective GV critical points, 
(b) Normalized consensus time as a function of x for the same 
models. 



drift is given by 
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where h is the local field (assuming even values be- 
tween -4 and 4), Prob(-l,/i) = (1 - x)Qx k (l - x) 4 - k , 
Prob(+l, h) = xQx k (l - x) 4 ~ k , and k = (4 + h)/2. In 
the case of the NLV model defined in Ref. [TT] , the flip- 
ping probability of a +1 spin takes the form /+(1) = r_4, 
/+(3/4) = r_ 2 , /+(l/2) = r , /+(l/4) = r 2 , and 
/+(0) = 0. Performing the summations in Eq. ^ we 
obtain the drift as a function of magnetization 



VM0)) t=o = ^(l-0 2 )F(0), 



(6) 



where 
F(</>) = 2(l-0 2 )(2r_ 2 



3r o ) + (3 + 2 )(r_ 4 -4r 2 ). (7) 



Depending on whether F((f>) is positive or negative, the 
drift will have the same sign of <fi or the opposite. For 
the NLV model, we have r 2 = r_4/4; therefore, in this 
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FIG. 3. (Color online) Average magnetization <j>(t) as a func- 
tion of time for different system sizes in the NLV model and 
the q — 4 qV model at the respective GV critical points, 
for initial conditions with magnetization 0(0) = —0.5, corre- 
sponding to x — 0.25. 
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FIG. 4. (Color online) Average drift as a function of the mag- 
netization for different system sizes in the NLV model in d = 2 
at the critical point, for system starting at an initial magne- 
tization 0(0) = -0.5 {x = 0.25). Top: t < 100. Bottom: 
t > 100. 



case F(4>) = 2(1 — </> 2 )(2r_ 2 — 3r ) and it is negative for 
r_ 2 < 3ro/2. For the values chosen for the NLV model, 
this inequality is indeed satisfied, and therefore the initial 
drift is positive for <fi < (x < 1/2), negative for <fr > 
(x > 1/2) and vanishes for <j> = (x = 1/2). 

These results clarify the origin of the nonlinear exit 
probability E(x) in .^-symmetric models. Initial con- 
ditions x 7^ 1/2 imply a strong nonzero drift, which 
rapidly brings the fraction of initial spins from its ini- 
tial value i to a different value x' . After this short 
transient the build up of spatial correlations cancels the 
drift, magnetization is conserved and the dynamics be- 
comes identical to that of linear VM. As a consequence 
Enlv(x) = Evm(x') = x' as witnessed by Fig. [3] where 
the density of +1 spins, starting from initial conditions 
x = 0.25, reaches a plateau x' = (1 + 4' 00 )/2 ~ 0.315, in 
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good agreement with the estimate of the exit probability, 
i.e. £(0.25) ~ 0.32. A similar mechanism is at the ori- 
gin of nonlinear exit probabilities for opinion dynamics 
models in d = 1 [25]- The same argument allows also 
to estimate the deviation of the consensus time from the 
entropic form Eq. ([!]). In the short initial transient the 
fraction of +1 spins quickly converges to x' — E(x). The 
consensus time is essentially set by the subsequent slow 
ordering, which occurs as in the linear VM, hence we can 
write 

T NLV (x)=T VM [E(x)}. (8) 

Figure [2]d) confirms the correctness of this theoretical 
estimate, and hints towards the relevance of the micro- 
scopic details of the model, which induce a strong drift at 
short time scales and lead in this way to nonuniversal fea- 
tures such as a nonlinear exit probability and anomalous 
consensus time. 



III. THE (j- VOTER MODEL 

In order to confirm the departure of the exit probabil- 
ity and of the consensus time from linear voter behavior 
for ^-symmetric models in the GV class, we consider the 
recently introduced non-linear g-voter (qV) spin model, 
which is defined as follows [2JJ: At each time step t, 
a random site i is chosen; additionally q nearest neigh- 
bor sites of i are also randomly selected (allowing for 
repetition to simplify the analysis), and their spins ex- 
amined. If all the q neighbors are in the same state, 
the spin at site i takes their common value. Otherwise, 
the spin at i flips its state with probability s. In any 
case, time is updated t — >• t + 1/N. With this defini- 
tion the flipping probability in the qV model takes the 
form f(xi,q,s) = x\ + e [1 — x\ — (1 — Xi) q ], where we 
remind that Xi is the fraction of neighbors of site i in the 
opposite state. Following simple mean-field arguments 
[TBI [T71 [2Tj , one can show the existence of a critical point 
e = e c (q), corresponding to GV behavior, separating a 
paramagnetic (disordered) phase for e > e c from a ferro- 
magnetic (ordered) phase at e < e c . In a d = 1 lattice, 
the qV model can be exactly mapped to the model of 
nonconservative voters proposed in Ref. |26j . From here, 
one observe voter behavior at e c — 1/2 for any value of q, 
while values of e ^ 1/2 lead to ordering dynamics with 
a non-trivial, non-linear exit probability. In the more 
interesting case d — 2, numerical evidence presented in 
Ref. [UJ for the case q = 4 indicated the presence of 
a critical point at e c ~ 1/4. For this value of e, evi- 
dence of voter behavior was found in terms of the decay 
of the density of interfaces and the scaling of the corre- 
lation function, both of which are fully compatible with 
the VM results. 

Here we confirm the value of the critical point of q = 4 
qV model by performing a finite size scaling analysis of 
the exit probability at x = 0.25, following the lines of the 
analysis carried out for the NLV model. In Fig. [5] (main 



0.230 




io' io 2 



FIG. 5. (Color online) Main: Finite-size scaling of the exit 
probability at x — 0.25 as a function of system size in the 
q — 4 qV model in d — 2, computed over at least 5 x 10 5 
independent realizations of the dynamics for each system size. 
Inset: Susceptibility and correlation length as a function of 
A = \e c — e\ for the q = 4 qV model on lattices of size L — 100 
and L = 5000, respectively. 

plot) we report the exit probability for this value of x and 
different values of e as a function of the lattice size L. A 
plateau is obtained for e ~ 0.24985, while larger (smaller) 
values of e lead to an increase (decrease) of E(x = 0.25) 
with L. We conclude that the critical point of the GV 
transition in d = 2 for the q = 4 qV model is located 
at e c — 0.24985(5), in good agreement with the previous 
estimate in Ref. [5TJ. We further confirm the fact that 
the critical q = 4 qV model belongs to the GV class by 
computing the exponents v and 7, Fig. [5] (inset). We 
obtain the values v ~ 0.57 and 7 ~ 1.28, in reasonable 
agreement with our estimates for the NLV model. 

As in the case of the NLV model, the exit probability 
at x = 0.25 indicates the presence of a nonlinear form, 
taking a value E(x = 0.25) = 0.222(7), smaller than 0.25 
beyond the estimated error bars. The non-linearity of 
E{x) is further confirmed in Fig. [2^,), where we compare 
the function E{x) for x < 0.5 with the linear form valid 
for linear VM. E{x) is non-linear in the whole range of 
x values, being independent of L at this critical point e c . 
The deviation of the qV model from linear VM behavior 
extends, similarly to the NLV model, to the functional 
dependence with x of the consensus time T/v(x) at the 
critical point e c , as shown in Fig.[2j}). 

Noticeably, in the case of the qV model, the exit prob- 
ability is smaller than x, in opposition to the NLV model, 
where we observed values E(x) > x. This smaller value 
of the exit probability is reflected in the evolution of the 
magnetization 0(f), see Fig. |3j which is again strongly 
not conserved at short times, exhibiting a sharp drop 
until it stabilizes, for times t > 50, at a plateau with 
approximate value (p^ rs —0.55. The time evolution of 
the average magnetization is again related with the drift 
v(4>). In Fig. [6] we plot (v(<j))) vs (f> for different system 
sizes in two distinct temporal regimes. For short times 
(t < 50) a sharp dip is present in the vicinity of the ini- 
tial magnetization. This is responsible for the initial de- 
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FIG. 6. (Color online) Average drift as a function of the 
magnetization for different system sizes in the q — 4 voter 
model in d — 2 at the critical point, for system starting at an 
initial magnetization <f>(0) = —0.5 (x — 0.25). Top: t < 50. 
Bottom: t > 50. 



crease of magnetization until it reaches the steady state 
and its conserved (in average) value. For larger values 
of time (t > 50) this dip is absent, and the drift takes 
a flat, almost vanishing form, thus ensuring conservation 
of magnetization. 

The dip in the drift a short times can also be under- 
stood by computing the drift in the initial uncorrelated 
condition. From Eqs. ^ and Q, and considering that 
for the qV model the flipping rates r, can be written as 
r_ 4 = 1, r_ 2 = 3(27 + 58e)/256, r = (1 + 14e)/16, 
r 2 = (1 + 174e)/256, and r± = 1, we obtain 

f(0) = ^[(l-2e)0 2 +41-114 £ ], (9) 

which is a positive function for all <j> for e < 41/114. 
Thus, for e — e c , we find that the initial drift is negative 
for 4> < (x < 1/2), positive for (j> > (x > 1/2) and 
vanishes for <f> = (x = 1/2). Again, for the qV model 
the nonlinear exit probability and the anomalous consen- 
sus time can be related through the argument leading to 
Eq. Q. Indeed, for x = 0.25, from Fig. [3] we read a 
stationary large time magnetization ip^ s» —0.55, corre- 
sponding to x' ~ 0.225, is in good agreement with the 
estimate of the exit probability, i.e. £7(0.25) ~ 0.222. 
Eq. ^ is again valid for the whole range of values of x, 
as shown in Fig. [2j 
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FIG. 7. (Color online) Main plot: Consensus time as a func- 
tion of x for different system sizes in the KKE model at 
the derealization transition in d — 2. Inset: Evolution of 
Tn(x) /Tn (0.5) for x = 0.25 as a function of the system size 
L. 



exit probability and a "non-entropic" form of the con- 
sensus time Tjv(x). By the same token, it is clear that 
in other type of models belonging to the GV class, con- 
servation of the average magnetization will guarantee a 
linear exit probability [4], as for linear VM. The question 
naturally arises about the form of the consensus times in 
those models in the GV class in which conservation of 
magnetization is enforced. 

To answer this question we have considered the KKE 
interface model at the delocalization transition |22j . 
which can be formulated in terms of a spin systems as fol- 
lows: At each time step a spin Sjj and one of its neighbors 
are randomly selected; if they are equal nothing happens, 
otherwise the number n + of positive neighbors of the neg- 
ative spin is computed and, with probability l/n + , Sij 
and the neighbor are made equal. This model belongs 
to the GV class, as has been shown in Ref. [IT], where 
a logarithmic decay of the density of interfaces was ob- 
served. In Fig. [7] we plot the rescaled consensus time as 
a function of x for different system sizes. These results 
indicate that again the consensus time deviates from the 
form expected in linear VM. The deviation from the be- 
havior given by Eq. should be attributed again to 
relevant microscopic details of the model. However, the 
mechanism inducing the deviation is here necessarily dif- 
ferent from the one responsible for the deviation for the 
NLV and qV models (i.e. initial short timescale noncon- 
servation of magnetization) . Its investigation constitutes 
an interesting line for future research. 



IV. CONSERVED MAGNETIZATION: THE 
KAYA, KABAKCIOGLU, AND ERZAN MODEL 



CONCLUSIONS 



From the analysis of the results obtained for the NLV 
and qV models, we expect that any ^-symmetric model 
belonging to the GV class and endowed with a nonlinear 
flipping probability f(xi) will exhibit a strictly nonlinear 



In summary, in this paper we have shown, using large 
scale numerical simulations of three different spin mod- 
els, that different subclasses of the generalized voter class 
for ordering dynamics actually exhibit different behav- 
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iors when unbalanced initial conditions are considered. 
All elements of the GV class are broadly characterized 
by lack of surface tension and a logarithmic decay of the 
density of interfaces in d = 2. However, when looking at 
the exit probability and the consensus time, two types of 
behavior occur. In one subclass, encompassing systems, 
such as KKE, with no Z 2 symmetry but magnetization 
conserved on average, the consensus time differs from 
the entropic form characterizing the VM, while the exit 
probability is linear. In the other subclass, composed 
by systems with Z 2 symmetry and nonlinear form of the 
flipping probability f(x), both the exit probability and 
the consensus time differ from the VM. This variation 
is essentially due to a nonconservation of magnetization 
during a short initial transient. The buildup of spatial 
correlations rapidly leads to an effective cancellation of 
the drift, so that the subsequent evolution is the same 
as for the linear VM, but starting from x 1 ^ x so that 
E(x) and T(x) are modified. To the best of our knowl- 
edge the results for the NLV and qV models reported 
here constitute the first example of a dynamics in d = 2 
with an exit probability different (in the large system size 



limit) from the step-function (typical of dynamics driven 
by surface tension) or the linear shape of VM. Nontrivial 
shapes were previously found only in d = 1 and 
in mean-field, but not in d = 2. Finally, our numerical 
estimates of the critical exponents associated to the GV 
transition suggest a value of the 7 exponent possibly dif- 
ferent from the mean-field value previously supposed to 
hold. Further research work is needed to fully clarify this 
issue. 
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